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Abstract 

I show that there exist universal constants C(r) < oo such that, for all 
loopless graphs G of maximum degree < r, the zeros (real or complex) of 
the chromatic polynomial Pg(q) he in the disc \q\ < C(r). Furthermore, 
C{r) < 7.963907r. This result is a corollary of a more general result on the 
zeros of the Potts-model partition function Zc(q, {v e }) in the complex anti- 
ferromagnetic regime |1 + v e \ < 1. The proof is based on a transformation of 
the Whitney-Tutte-Fortuin-Kasteleyn representation of Zc(q, {v e }) to a poly- 
mer gas, followed by verification of the Dobrushin-Kotecky-Preiss condition 
for nonvanishing of a polymer-model partition function. I also show that, for 
all loopless graphs G of second-largest degree < r, the zeros of Pg(q) he in the 
disc \q\ < C(r) + 1. Along the way, I give a simple proof of a generalized (mul- 
tivariate) Brown-Colbourn conjecture on the zeros of the reliability polynomial 
for the special case of series-parallel graphs. 

KEY WORDS: Graph, maximum degree, second-largest degree, chromatic polyno- 
mial, dichromatic polynomial, Whitney rank function, Tutte polynomial, reliability 
polynomial, chromatic roots, Potts model, antiferromagnetic Potts model, Fortuin- 
Kasteleyn representation, polymer gas, Mayer expansion, cluster expansion, phase 
transition. 



1 Introduction 



The polynomials studied in this paper arise independently in graph theory and 
in statistical mechanics. It is appropriate, therefore, to begin by explaining each of 
these contexts. Specialists in these fields are warned that they will find at least one 
(and perhaps both) of these summaries excruciatingly boring; they can skip them. 

Let G = (V, E) be a finite undirected graph[| with vertex set V and edge set E. 
For each positive integer q, let Pg{<i) be the number of ways that the vertices of G can 
be assigned "colors" from the set {1,2, . . . , q} in such a way that adjacent vertices 
always receive different colors. It is not hard to show (see below) that Pg{q) is the 
restriction to Z + of a polynomial in q. This (obviously unique) polynomial is called 
the chromatic polynomial of G, and can be taken as the definition of Pg{q) for 
arbitrary real or complex values of q.f\ 

The chromatic polynomial was introduced in 1912 by Birkhoff |13|]. The original 



hope was that study of the real or complex zeros of Pg(q) might lead to an analytic 
proof of the Four-Color Conjecture j75], j38|l , which states that -Pg(4) > for all loopless 
planar graphs G. To date this hope has not been realized, although combinatoric 
proofs of the Four-Color Theorem have been found |2|, |3], f|, |112|| . Even so, the 
zeros of Pg{(i) are interesting in their own right and have been extensively studied. 
Most of the available theorems concern real zeros g^, [TT7|, [128], |12|, [56], [130], |TT3|, [38 



but there has been some study (mostly numerical) of complex zeros as well [ 52J, [10 



A more general polynomial can be obtained as follows: Assign to each edge e G E 
a real or complex weight v e . Then define 

Z G (q, R}) = II l 1 + V e$( (r x 1 (e),C- X2ie )) , (1.1) 

{a x } eeE 

where the sum runs over all maps a: V — > {1, 2, . . . , q}, the 5 is the Kronecker delta, 
and Xi(e) , X2(e) G V are the two endpoints of the edge e (in arbitrary order). It is not 
hard to show (see below) that Zc(q, {v e }) is the restriction to q G Z + of a polynomial 
in q and {v e }. If we take v e = —1 for all e, this reduces to the chromatic polynomial. 
If we take v e = v for all e, this defines a two- variable polynomial Zc{q,v) that was 
introduced implicitly by Whitney ||124| , |125| , |126j| and explicitly by Tutte [|115| , |116|| ; 



it is known variously (modulo trivial changes of variable) as the dichromatic poly- 
nomial, the dichromate, the Whitney rank function or the Tutte polynomial 

1 In this paper a "graph" is allowed to have loops and/or multiple edges unless explicitly stated 
otherwise. 
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Two excellent reviews on chromatic polynomials are [plL B3J. An extensive bibliography on 



chromatic polynomials is 1 30 
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GUM 

In statistical mechanics, (|1 . 1| ) is known as the partition function of the qr-state 
Potts model.Q In the Potts model |13T], |132[, an "atom" (or "spin") at site x E V 



can exist in any one of q different states (where q is an integer > 1). The energy 
of a configuration is the sum, over all edges e £ E, of if the spins at the two 
endpoints of that edge are unequal and — J e if they are equal. The Boltzmann 
weight of a configuration is then e~^ H , where H is the energy of the configuration 
and (3 > is the inverse temperature. The partition function is the sum, over all 
configurations, of their Boltzmann weights. Clearly this is just a rephrasing of (|1.1|) , 
with v e = e@ Je — 1. A coupling J e (or v e ) is called ferromagnetic if J e > {y e > 0) 
and antiferromagnetic if — oo < J e < (— 1 < v e < 0). 

To see that Zc{q, {v e }) is indeed a polynomial in its arguments (with coefficients 
that are in fact or 1), we proceed as follows: In (|1.1|) , expand out the product 
over e 6 £, and let E' C E be the set of edges for which the term v e 5 a , , a , , is 
taken. Now perform the sum over configurations {a x }: in each connected component 
of the subgraph (V, E') the spin value a x must be constant, and there are no other 
constraints. Therefore, 

Z G (q,{v e }) = J> W n^' (1-2) 

E'CE e£E' 

where k(E') is the number of connected components (including isolated vertices) 



in the subgraph (V,E'). The expansion (1^2) was discovered by Birkhoff ]T3[ and 



Whitney [ |124j| for the special case v e = — 1 (see also Tutte ||115| , |116|| ); in its general 
form it is due to Fortuin and Kasteleyn 0, (see also P9| ). We take ( |1.2| ) as the 
definition of Zc{q, {v e }) for arbitrary complex q and {v e }. 

In statistical mechanics, a very important role is played by the complex zeros of 
the partition function. This arises as follows [|133|| : Statistical physicists are inter- 
ested in phase transitions, namely in points where one or more physical quantities 
(e.g. the energy or the magnetization) depend nonanalytically (in many cases even 
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The Tutte polynomial Ta(x,y) is conventionally denned as [123, p. 45] flllL pp. 73, 101] 



T G (x,y) = ^(aj-iJ^-M^^-ijlB'l+^'Hvl 



E'CE 



where k(E') is the number of connected components in the subgraph (V, E'). Comparison with (1.2) 
below yields 

T G (x,y) = (x-l)- k W(y-l)-Wz G ((x-l)(y-l),y-l) . 



The Potts model |80 



was invented in the early 1950's by Domb (see 
5l, was invented in 1920 by Lenz j7Cj (see |2l|, |63| 



The q = 2 case, 
091). The q = 4 case, 



known as the Ising model 

which is a special case of the Ashkin-Teller model, was invented in 1943 by Ashkin and Teller j|] 
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discontinuously) on one or more control parameters (e.g. the temperature or the mag- 
netic field) . Now, such nonanalyticity is manifestly impossible in (|1.1|) / (|1.2p for any 
finite graph G. Rather, phase transitions arise only in the infinite-volume limit. That 
is, we consider some countably infinite graph G^ = (V^^oo) — usually a regular 
lattice, such as Z d with nearest-neighbor edges — and an increasing sequence of finite 
subgraphs G n = (V n ,E n ). It can then be shown (under modest hypotheses on the 
G n ) that the (limiting) free energy per unit volume 

f Goo (q,v) = lim \Vn\~ 1 log Z Gn (q,v) (1.3) 

n— >oo 

exists for all nondegenerate physical values of the parameters^], namely either 

(a) q integer > 1 and — 1 < v < oo [using (|1.1|) : see e.g. |55], Section 
1.2]] 

or (b) q real > and < v < oo [using (|1.2|) : see 

This limit fG ao ( ( l, v ) is m general a continuous function of v; but it can fail to be a 
real-analytic function of v, because complex singularities of log Z Gn (q, v) — namely, 
complex zeros of Zc„{q, v) — can approach the real axis in the limit n — > oo. There- 
fore, the possible points of physical phase transitions are precisely the real limit points 
of such complex zeros. As a result, theorems that constrain the possible location of 
complex zeros of the partition function are of great interest. In particular, theo- 
rems guaranteeing that a certain complex domain is free of zeros are often known as 
Lee-Yang theorems^ 

The purpose of this paper is to prove an upper bound on the complex g-plane 
zeros of the Potts-model partition function Zc{q, {v e }), valid throughout the "com- 
plex antiferromagnetic regime" |1 + v e \ < 1, under certain "local" conditions on the 
weights {v e }: for example, in terms of the quantity max ^2 \v e \. As a corollary, I 

obtain upper bounds on the zeros of the chromatic polynomial Pg{q) m terms of the 
maximum degree of the graph G. More precisely, I show that there exist universal 
constants C(r) < oo such that, for all loopless graphs G of maximum degree < r, 
the zeros of Pg{q) he in the disc |g| < C(r). This answers in the affirmative a ques- 
tion posed by Brenti, Royle and Wagner Question 6.1], generalizing an earlier 
conjecture of Biggs, Damerell and Sands |T2| limited to r-regular graphs. The con- 
stants C(r) arise as the solution of an explicit minimization problem, and I prove 

5 Here "physical" means that the weights are nonnegative, so that the model has a probabilistic 
interpretation; and "nondegenerate" means that we exclude the limiting cases v = — 1 in (a) and 
q = in (b), which cause difficulties due to the existence of configurations having zero weight. 

6 The first such theorem, concerning the behavior of the ferromagnetic Ising model at complex 
magnetic field, was proven by Lee and Yang |6!| in 1952. A partial bibliography (through 1980) of 
generalizations of this result can be found in f[M. 



51, Theorem 4.1] and 
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that C(r) < 7.963907r. This linear dependence on r is best possible, as the example 
of the complete graph K r+ i shows that C(r) > r. 

Furthermore, I show that the presence of one vertex of large degree cannot lead to 
large chromatic roots. More precisely, if all but one of the vertices of G have degree 
< r, then the zeros of Pg(q) he in the disc |g| < C{r) + 1. Please note that a result 
of this kind cannot hold if "all but one" is replaced by "all but two" , for in this case 
the chromatic roots can be unbounded, even when r = 2 and G is planar [103 . 



The proofs of these results are based on well-known methods of mathematical 
statistical mechanics. The first step is to transform the Whitney-Tutte-Fortuin- 
Kasteleyn representation (|1.2|) into a gas of "polymers" interacting via a hard-core 
exclusion (Section^). I then invoke the Dobrushin condition |34], (or the closely 
related Kotecky-Preiss condition [^, |104|| ) for the nonvanishing of a polymer- model 
partition function (Section ^). Lastly, I verify these conditions for our particular 
polymer model, using a series of simple combinatorial lemmas, some of which may be 
of independent interest (Section in particular, I give a simple proof of a generalized 
(multivariate) Brown-Colbourn conjecture on the zeros of the reliability polynomial 
for the special case of series-parallel graphs (Remark 3 in Section |4.1|) . The main 
results of this paper are contained in Section some generalizations and extensions 
are in Section ||. I conclude with some conjectures and open questions (Section |7|). 

With a little more work, it should be possible to extend the arguments of this 
paper to prove the existence and analyticity of the limiting free energy per unit volume 
( |1.3|) for suitable regular lattices G^ and translation-invariant edge weights v e , in the 
same region of complex q- and {t> e }-space where Z will be proven (in Section [5]) to 
be nonvanishing uniformly in the finite subgraphs V n ("uniformly in the volume" in 
statistical-mechanical lingo). In particular, this would provide a convergent expansion 
for the limiting free energy in powers of 1/q. However, I have not worked out the 
details. 

This paper would never have seen the light of day without the help and advice 
of Antti Kupiainen. During my visit to Helsinki in September-October 1997, I told 
Antti of my conjectures about -Pg(q') an d Z G (q, {v e }) — conjectures that I had no 
good idea how to prove. He immediately saw that they ought to be provable by 
cluster (or Mayer) expansion. My reaction was, "Ugh! You know how I detest the 
cluster expansion!"; indeed, I had resisted learning it for nearly 20 years and had 
devoted much of my work in mathematical physics to finding ways of circumventing 
it [ 102 , |23| , [24] , Antti assured me that the cluster expansion is not so difficult, and 
he suggested that I study the excellent review article of Brydges |2~2] . We also quickly 



figured out how to represent Zc(q, {ve}) as a polymer gas. Jean Bricmont then told 
me about the work of Kotecky and Preiss f65 |, and Roman Kotecky informed me of 
the work of Dobrushin [[53] . Here, finally, was a version of the cluster expansion simple 
enough that even I could understand it! Nine months later, I figured out how to verify 
the Dobrushin (or Kotecky-Preiss) condition and thereby complete the proof. 
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2 Transformation of the Potts- Model Partition Func- 
tion to a Polymer Gas 

Let G = (V, E) be a finite undirected graph equipped with complex edge weights 
{v e } e< zE- If G contains a loop e (i.e. an edge connecting a vertex to itself), this simply 
multiplies Z G (q, {v e }) by a factor 1 + v e ; so we can assume without loss of generality 
that G is loopless, and we shall do so in this section in order to avoid unnecessary 
complications. Likewise, if G contains multiple edges ei, . . . , e n connecting the same 
pair of vertices, they can be replaced, without changing the value of Z, by a single 
edge e with weight v e = Yli=i(^ + u eJ — 1- So we could assume without loss of 
generality, if we wanted, that G has no multiple edges. But this assumption would not 
simplify most of our subsequent arguments, so we shall usually refrain from making 
it. Note, however, that our numerical bounds frequently get better if multiple edges 
are replaced by a single equivalent edge. 

So let G be loopless, and consider the Whitney-Tutte-Fortuin-Kasteleyn repre- 
sentation ( |1.2[ ) of the Potts-model partition function Z G {q, {v e }). For each term in 
( |1.2| ) we decompose the subgraph (V, E') into its connected components. Some of 
these components may consist of a single vertex and no edges; the remaining com- 
ponents are disjoint connected subgraphs (Si, Ex), . . . , (Sn,En) with \Si\ > 2. The 
total number of components is 

/ N \ 

(2.1a) 
(2.1b) 




It follows that: 



Proposition 2.1 (jointly with Antti Kupiainen) Let G = (V,E) be a loopless 
finite undirected graph equipped with edge weights {v e } e( zE- Then 



where 



and 



J polymer,G 



Zc(q, {v e }) = q lvl z polymertG (q, {v e }) , (2.2) 

oo N 
^=0 ' Si, ...,Sjv disjoint i=l 

q -(\s\-D E Uve if \S\>2 

if\S\<l 



. E C E e£E 

W [S) — \ (S, E) connected (2-4) 
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The sum in (2J5) runs over pairwise disjoint subsets S±, . . . , Sn of V , and the term 



N = in ( \2.3j ) is understood to contribute 1. 



Note, in particular, that w(S) — if S is disconnected [i.e. if the induced subgraph 
(S,Es) is disconnected]. 



The "polymer model" (|2.3|) - (|2.4j) has the form of a grand-canonical gas (see Sec- 
tion |3] for the precise definition) 

oo j N 

Z polymer , G (q,{v e }) = U w ^ II W( S t> S i) ( 2 - 5 ) 

N=0 ' S!,...,S N i=l l<i<j<N 

with single-particle state space V*(V) [the set of all nonempty subsets of V), fugacities 
w(S), and two-particle Boltzmann factor given by a hard-core exclusion 

W(S,S') = ( 1 if SnS ' = ® (2.6) 
I otherwise 



Graph theorists will recognize the right-hand side of ([2.5) as the generating function, 



in the variables w(S), for independent subsets of vertices of the intersection graph of 
V*{V). 

The usefulness of (|2.2|)-(|2T6|) comes from the fact that the fugacities ^(5*) are all 
suppressed by powers of hence are small for large \q\. Moreover, if the sum over 
E in ([2.4|) can be controlled, one expects that w(S) will be exponentially decaying 
in \S\ when \q\ is large enough. This raises the hope that the Mayer expansion [|119|[ , 



which is an expansion of log Z po i ymer! G in powers of the fugacities u'(S'), might converge 
for sufficiently large \q\. If so, this would imply that Z po iy merG ^ in the region of 
convergence. That is what we go about proving in the following sections — but in 
the opposite order. 



3 Dobrushin and Kotecky— Preiss Conditions for 
the Nonvanishing of Z 

In statistical mechanics, a grand-canonical gas is defined by a single-particle 
state space X (here assumed for simplicity to be finite), a fugacity vector w = 
{w x } x& x £ C x , and a two-particle Boltzmann factor W(x,y) [a symmetric func- 
tion W: X x X — > C]. The (grand) partition function Z(w,W) is then defined 
to be the sum over ways of placing N > "particles" on "sites" xi,...,xjsr G X, 
with each configuration assigned a "Boltzmann weight" given by the product of the 
corresponding factors w Xi and W(xi,Xj): 

oo j N 

z(w,w) = j2m s n wfaxj), (3.i) 

N=0 ' x 1 ,...,x N £X i=l l<i<j<N 
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where the N = term is understood to contribute 1. Under very mild conditions 
on W [e.g. |W(a;, y)\ < 1 for all x, y is more than sufficient], Z(w,W) is an entire 
analytic function of w. Our goal is to find a sufficient condition for Z(w, W) to be 
nonvanishing in a polydisc Dr = {w: \w x \ < R x }. This would imply, in particular, 
that log Z(w, W) is an analytic function of w in Dr. 
We say that W is 

• physical if < W(x, y) < +oo for all x, y G X 

• repulsive if |W(a;, < 1 for all x,y G X 

• physical and repulsive if < W(x, y) < 1 for all x, y G X 

• hard-core if W(x, y) = or 1 for all x, y G X 

• hard-core self -repulsive if VK(x, x) = for all x G X 

An important special case is when W is hard-core and hard-core self-repulsive: then 
Z(w, W) is the generating function for independent sets of vertices of the graph 
G = (X, E) defined by placing an edge between each pair of vertices x ^ y for which 
W{x,y) = 0. 



Dobrushin |53, has given an elegant sufficient condition for the nonvanishing 
of Z in a polydisc Dr, whenever W is hard-core and hard-core self- repulsive. His 
proof is astoundingly simple, avoiding all the combinatoric complication that has 
given cluster expansions such a reputation for difficulty. Here I shall present a slight 
extension of Dobrushin's theorem, in which the condition of hard-core interaction 
is replaced by the weaker assumption that the interaction is physical and repulsive; 
moreover, the conclusion of the theorem is slightly strengthened. (We won't really 
need this extension — the original Dobrushin theorem would suffice for our purposes 
- but the stronger result is no more difficult, and it gives a bit more insight into 
the method of proof.) The hard-core se//-repulsion is, however, essential both in 
Dobrushin's version and in my own: it guarantees that each "site" x G X can be 
occupied by at most one "particle" It follows that the partition function can be 
rewritten as a sum over subsets: 

z( W ,w) = j2 n>* n w ( x >y)> ( 3 - 2 ) 

X'CX x£X> (xy)eX' 

where the second product runs over unordered pairs x, y G X' [x ^ y) with each pair 
counted once. 

Let us define, for each subset A C X, the restricted partition function 

z A ( W ,w) = j2 n>* n w &v)- ( 3 - 3 ) 

X'CA xeX' (xy)eX' 



s 



Of course this notation is redundant, since the same effect can be obtained by setting 
w x = for x G X \ A, but it is useful for the purposes of the inductive proof. We 
have: 



Theorem 3.1 Let X be a finite set, and let W satisfy 

(a) < W{x, y) < 1 for all x, y G X 

(b) W(x,x) = for allxeX 

Suppose there exist constants R x > and < K x < 1/R X satisfying 

l-W(x,y)K y R y 



K x > n 

y+x 



1 - KyRy 



(3.4) 



for all x G X . Then, for each subset ACI, Z A (w, W) is nonvanishing in the closed 
polydisc D R = {w G C x : \w x \ < R x } and satisfies there 



d log Z A (w, W) 



K x 



<9uu 



for all x £ A 



< J 1 - K x \w x \ 

[o for allx G X \ A 

Moreover, if w, w' G D R and w'Jw x G [0, +oo] for each x G A, then 

Za(V,W) 



(3.5) 



log 



Z A (w,W) 



<E 

xeA 



log 



1-K x \w' x \ 



1 — K x \w x \ 



where on the left-hand side we take the standard branch of the log, i.e. \ Imlog • 



(3.6) 



< 7T. 



Remarks. 1. It follows from (|3.4j) that K x > 1 and hence that R x < 1. 



2. The conclusion of Dobrushin's theorem |33, |3^] is the special case of ( |3.6| ) in 
which some of the w' x are equal to w x and others are equal to 0, and in which only 
the real part of the logarithm on the left-hand side is handled. 



Proof. Note first that (|3.5|) for any given A implies (3J;) for the same A, by inte- 
gration. 

The proof is by induction on the cardinality of A. If A = the claims are trivial. 
So let us assume that ( |3.5|) [and hence also (|3.6|) 1 holds for all sets of cardinality 
< n, and let a set A of cardinality n be given. Let x be any element of A, and let 
A' = A \ {x}. It follows from flO) that 



Z A (w,W) = Z A ,(w,W) + w x Z Af (w,W) 



(3.7) 
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where 

w y = W(x,y)w y ; (3.8) 
here the first term on the right-hand side of (|3.7|) covers the summands X' ^ x, while 



the second covers X' 3 x. Note that w G Dr since |W(x, y)\ < 1. From ( |3.7| ) we have 

k{w) 



d 

dw x 



\ogZ K {w,W) 



where 



k{w) 



1 + k(w)w x 
Z A >(w,W) 



(3.9) 



Now by the inductive hypothesis 
W(x, y) > 0, we have 



lHw)l < JJ i-^(x,y)^K 



Z A >(w,W) ' 
for A', and using the fact that w y /w y 



(3.10) 



yeA' 



y\ UJ y\ 



s n 

yex\{x} 



l-W{x,y)K y 


Kl 


1-Ky 


Wy\ 





(3.11) 



which is < K x by the hypothesis ( |3.4j ). This proves ( |3.5| ) for A, and hence completes 
the induction. | 

Let us now return to the special case of a hard-core interaction. If W(x, y) = 
(resp. 1), we say that x and y are incompatible (resp. compatible) and write x ^ y 
(resp. x ~ y). Note that in our convention x ^ x, in agreement with some authors' 



convention |65|, |9l], [L01|| and contrary to others' |34| , p5[ . The hypothesis ( p.4|) is then 
equivalent to the existence of constants c x > such that 



Rr. < (e c 



1) exp i - c y j 



(3.12) 



for all x G X [set c x = — log(l — K X R X )]. This is the Dobrushin |33|, ^] condi- 
tion. Slightly stronger, and more convenient to check, is the Kotecky-Preiss [^, |104 | 
condition 



R x < c x exp I - ^ c y | 



(3.13) 



Let us now consider the important special case in which the single-particle state 

oo 

space X can be partitioned as X = |J X n in such a way that 



n=l 



i? y < v4„m for all x G X„ 

y&X n : y^x 



(3.14) 
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for suitable constants {^4n}£!Li- [This typically arises when X is some set of nonempty 
subsets of a finite set V, and x ^ y means x fl y ^ 0; we will then take X n to be the 
sets of cardinality n, and will prove ( |3.14| ) by proving 



R y < A n for alH £ V , (3.15) 



which is manifestly stronger than ( 3.14)) .] Let us take 



c x = e an R x for all x £ X n (3.16) 
with some suitably chosen a > 0. Then, for ( |3.14j ) to imply the Kotecky-Preiss 



condition (|3.13|) , it suffices that 

oo 

J2e an A n < a. (3.17) 



n=l 



We have therefore proven: 

oo 

Proposition 3.2 Suppose that X = [J X n (disjoint union) and that there exist 
constants {An} 1 ^^ and a such that 



n=l 



(a) Ry — A n m for all m,n and all x £ X„ 

yeX n : yfix 
oo 

(b) J2e an A n < a. 

n=l 



Then the Kotecky-Preiss condition ( S.lSj) holds with the choice c x = e an R x for x £ 



Remarks. 1. Suppose we try the more general Ansatz c x = b n R x for x £ X n . 
Then ( |3.14| ) implies the Kotecky-Preiss condition ( 3.13| ) in case b n > e an where 
a = Y^nLi b n A n . But in that case b' n = e an > e an where a' = Yl^Li b' n A n . So there 
is no loss of generality in restricting attention to b n = e an for some a. 

2. Since the state space X is finite, only finitely many of the A n are nonzero. 
Nevertheless, we often have occasion to consider simultaneously an infinite family of 
problems — e.g. in this paper, all loopless graphs G of maximum degree < r and 
arbitrarily many vertices — and it is natural to seek bounds that are uniform over 
the family. So it is useful to forget that only finitely many of the A n are nonzero. 
(Moreover, similar methods can be applied to problems with an infinite state space 
X, in which case {A n } is a genuinely infinite sequence.) This leads to two further 
remarks: 
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3. For the condition 

oo 

3a > such that ^e an A n < a (3.18) 

n=l 

to hold, it is necessary that the sequence {Aj}^ have some exponential decay (i.e. 
A n < Ce~ en for some e > 0), but there is no minimum required rate of decay. Indeed, 
if {Ai}^Li has any exponential decay at all, then by modifying finitely many of the 
A n one can make ( |3.18| ) hold. It can thus be valuable in applications to work hard 
on estimating the first few coefficients A n (see |59|] for an example) .[] 

4. Let 5 = liminf n ^ 00 (— log A n )/n. Then F(a) = or 1 Y^=\ ean A n is finite-valued 
and continuous (in fact, real- analytic) on < a < 5, left- continuous (as a map into 
the extended real line) as a f 5, and identically +oo for a > 6. In particular, the 
infimum of F(a) is attained, so (|3.18| ) is equivalent to 

oo 

inf a~ 1 '^2e an A n < 1 . (3.19) 

n=l 



Important Final Remark. The results in this section provide an extraordinarily 
simple proof of the convergence of the Mayer expansion for a grand-canonical gas 
with physical and repulsive two-particle interactions. To see what is at issue, let us 
first trivially rewrite the partition function ( |3.1| ) as 

oo ^ N 

z(w,w) = j2m s ii^ e n f > ( 3 - 2 °) 

N=0 ' xi,...,x N £X i=l GeSiv {ij)eG 

where Qn is the set of all (simple loopless undirected) graphs on the vertex set 
{1, . .., N}, and 

F(x,y) = W(x,y)-1 (3.21) 

is called the two-particle Mayer factor. Then standard combinatorial arguments ||119|| 
show that 

oo - N 

log Z(w,W) = ^— II^. E II (3-22) 

JV=1 ' x lt ...,x N eX i=l GeC]v (ij)eG 

at least in the sense of formal power series in w, where Cn C Q n is the set of connected 
graphs on {1, ... , N}. This is the Mayer expansion; the principal problem is to prove 
its convergence in some specified polydisc. The usual approach to proving convergence 

7 The emphasis in S, |9l], [35| [lOlf on the special case A n = Ce~ en with C = 1 is thus 
somewhat misleading, inasmuch as it suggests that there is a minimum allowed rate of decay e. 
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of the Mayer expansion [79, 91, 28, 22 



101, 26, 104] is to explicitly bound 



Zg, [27|, 

the terms in Q3.22 ); this requires some rather nontrivial combinatorics (for example, 
Proposition ^J] below together with the counting of trees). Once this is done, an 
immediate consequence is that Z is nonvanishing in any polydisc where the series for 
logZ is convergent. Dobrushin's brilliant idea |34], |35| was to prove these two results 
in the opposite order. First one proves, by an elementary induction on the cardinality 
of the state space, that Z is nonvanishing in some specified polydisc (Theorem |3.1| ); 
it then follows immediately that log Z is analytic in that polydisc, and hence that its 
Taylor series (|3.22|) is convergent there. It is an interesting open question to know 
whether this approach can be made to work without the assumption of hard-core 
self-repulsion. 



4 Some Combinatorial Lemmas 



4.1 Reduction to trees 



The weight w(S) involves a sum ( |2.4|) over connected subgraphs (S, E) of the 
induced subgraph (S,Es). The trouble is that there may be "too many" connected 
subgraphs. It is remarkable, therefore, that this sum can sometimes be bounded by a 
sum over a much smaller set of graphs, namely spanning trees. The following propo- 
sition underlines the special role played by the "complex antiferromagnetic regime" 
A = {v E C: |l + v| < 1}. 

Proposition 4.1 (Penrose [ |79|| ) LetG= (V,E) be a finite undirected graph equipped 
with complex edge weights {v e } e€ E satisfying \l + v e \ < 1 for all e. Then 



e n 



E'CE 
(V, E') connected 



e£E' 



< 



e n 



(4.1) 



E' C E 
(V, E') tree 



edE' 



Penrose |79| proved this when G is the complete graph K n ; the result then follows for 
all graphs without loops or multiple edges (it suffices to set v e = on the nonexistent 
edges). Here I present a minor modification of Penrose's proof that permits loops and 
multiple edges: 

Proof. We can assume without loss of generality that G is connected, since otherwise 
both sides of the inequality are zero. Let C (resp. T) be the set of subsets E' C E 
such that (V,E') is connected (resp. is a tree). Clearly C is an increasing family 
of subsets of E with respect to set-theoretic inclusion, and the minimal elements 
of C are precisely those of T (i.e. the spanning trees). It is a nontrivial but well- 
known fact |IR Sections 7.2 and 7.3] ||134| , Section 8.3] that the (anti-) complex C is 
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partitionable: that is, there exists a map R: T — > C such that R(T) D T for all 
T G T and C = [+J [T, R(T)] (disjoint union), where [Ei,E 2 ] denotes the Boolean 
interval {E'\ E\ C £" C E 2 }. In fact, many alternative choices of R are available fll"5" 



Sections 7.2 and 7.3] j|7|, Sections 2 and 6] JTT|, Proposition 13.7 et seq.], and none 
of the subsequent arguments will depend on a specific choice of R. Nevertheless, for 
completeness, we shall give at the end of this proof a concrete construction of one 
possible R. 

Given the existence of R, we have the immediate identity 

e n«= = e n* e n * 

E' CE ee-E' T C E eeT TCE'CR(T) eeE'\T 

(V, E') connected (V,T) tree 

= e n** n (4.2) 

TCE e€R(T)\T 
(V,T) iree 



In particular, if |1 + v e \ < 1 for all e, then ( |4.1| ) follows. 

We now indicate a construction of R that is a slight variant of the one used by 
Penrose J70| (he orders the vertices, while I order the edges): Choose (arbitrarily) a 
vertex x G V and call it the root; and choose (arbitrarily) a numbering of the edges. 
For each E' G C and y G V, let depth E i(y) be the length of the shortest path in E' 
connecting y to the root. For each y G V \ {x}, let e(y) be the lowest-numbered 
edge in E' connecting y to a vertex y' with depth E ,(y') = depth E ,(y) — 1. And finally, 
let S(E') = {e(y): y G V \ {x} }. Then trivially S(-E') C moreover, it is easy 
to see that (V, S(E')) is a tree and that depth s ^ E ,^(y) = depth E ,(y) for all y G V. 
Conversely, given a spanning tree (V, T), it is not hard to see that S(E') = T if and 
only if T C E' C R(T), where R(T) is obtained from T by adjoining all edges e G E 
that 

(a) connect two vertices of equal depth T (this includes loops, if any), or 

(b) connect a vertex y to a vertex ?/' having depth T (y') = depthrp(y) — 1 where e 
is higher-numbered than the edge already in T that connects y to a vertex y" 
having depth T (y") = depth T (y) — 1. 

This completes the proof. ■ 

Remarks. 1. The identity ( f4.2| ) and the inequality ( [4.1| ) generalize to matroids. 
Indeed, for any matroid M, the independent sets of M form a simplicial complex 
IN(M), called a matroid complex; moreover, every matroid complex is shellable, and 



every shellable complex is partitionable [15, Theorem 7.3.3 and Proposition 7.2.2]. 



For the cographic matroid M*(G) } the independent sets are the complements of con- 
nected subgraphs, and the bases are complements of spanning trees, so we recover the 
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situation of Proposition O. I thank Criel Merino for teaching me about matroids 
and for helping me notice a silly error in my original proof of Proposition |4.1| . Earlier, 
Dave Wagner had informed me that analogues of ( }4.2| ) hold for shellable simplicial 
complexes (see [[107| , sections 0.3 and III. 2] for the definition). There is a long history 
of identities related to Q: see, for example, |87], [73], |22], |2|, [27|, 0, |26| and the 
references cited therein. 

2. I conjecture that ( f4.1p can be strengthened so that on the right-hand side the 
absolute value is put outside the sum rather than inside. (This would be useful in 
case the {v e } do not all have the same phase.) In fact, I conjecture more: Let 

c{E') = \E'\ - \V\ + k(E') (4.3) 

be the cyclomatic number of the subgraph (V,E r ), and define the generalized con- 
nected sum 

c G (x,{ve}) = Yi xcm a* ( 44a ) 

E' C E eeE' 
(V, E') connected 

= A 1 -l y IC G (l,{Ai; e }) (4.4b) 

= limg- 1 A- |y| Z G (Ag,{At; e }) . (4.4c) 

q->0 

In particular, A = corresponds to the tree sum and A = 1 to the connected sum. 
Then I conjecture that 

(a) If |1 + v e \ < 1 for all e, then |Cg(A, {f e })| is a decreasing function of 
A on < A < 1. 

I had originally conjectured a stronger result, namely 

(a;) If \r + v e \ < r for all e, then (-1)™ (d n /d\ n ) \C G (\, {v e })\ 2 > on 
< A < 1, for all n > 

either for r = 1 or, failing that, for r = |; but this is in fact false for all r > 0, even 
for the second derivative evaluated at A = with equal edge weights v e = v. Indeed, 
if we write 

C G (A, v) = v\ y \- 1 [oq + a x v\ + a 2 v 2 X 2 + . . .] (4.5) 
where aj is the number of spanning subgraphs of G having j cycles, then 



(d 2 /d\ 2 )\C G (\,v)f 



> holds for all v if a\ > 2aoQ2, but fails for v in a wedge 



A=0 

near the imaginary axis if a\ < 2a a 2 . Now the complete bipartite graph K 3 4 has 
a = 432, ai = 612, a 2 = 456 and hence provides a counterexample. Nevertheless, 
(a^.) might be true for some interesting subclasses of graphs G. 

For real v e G [—1,0], by contrast, I am able to prove a result even stronger than 
(a'^), namely 
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(b) If -1 < v e < for all e, then (-lf+l^l- 1 (d n /d\ n ) C G {\, {v e }) > 
on < A < 1, for all integers n > 0. 



Indeed, by fl4.2| ) and (4.4b), we have 
C (A,K}) = 



e n ( i+a ^) 

TCB eGT eeR(T)\T 
(V, T) free 



(4.6) 



and hence 
d 



d\ n 



E 



e n 



n 



(4.7) 



TP T C R(T) \ T eGTUT eeR(T)\(TUT) 
(V,T) free |f| = n 



which has the claimed sign whenever < A < 1 and — 1 < v e < for all e. 

3. C*g(1, {^e}) is equal, up to a prefactor, to the reliability polynomial Rc({Ve}) 



32| , where p e is the probability that edge e is operational and t> e = p e /(l — Ve)'- 



ee_B 



Ve 



C G (l,{ Pe /(l-Ve)}) 



(4- 



Now the Brown-Colbourn conjecture |2(], |120|| states that for any connected graph G 
(loops and multiple edges are allowed), Rg(v) whenever \p — 1| > 1. A more 
general conjecture is that Rc({Ve}) 7^ whenever \p e — 1| > 1 for all edges e, or 
equivalently, that (7^(1, {f e }) ^ whenever < |1 + v e \ < 1 for all e. But this 
generalized Brown-Colbourn conjecture is an immediate consequence of conjecture 
(a): for if we had Cq(1, {v e }) = with |1 + v e \ < 1 for all e, then we could choose 
e > such that v' e = (1 + e)v e satisfy |1 + v' e \ < 1 for all e, and we would have 
C G (X, {v' e }) = for A = l/(l+e) [but not, of course, identically for l/(l + e) < A < 1]. 

Note also that if the generalized Brown-Colbourn conjecture holds for a graph 
G, then it holds also for any graph that can be obtained from G by a sequence 
of doublings of edges ("parallel expansions") and/or subdivisions of edges ("series 
expansions"). This follows from the formulae |33, p. 35] 



RG'ttVe, Vl,V2}) = ifc({p e , l-(l- P i)(l-p 2 )}) 
R G >({Ve,Vl,V2}) = [l-(l-p 1 )(l-p 2 )] J R G ({p e , 



VlV2 



V1+V2- VlV2 



(4.9) 
(4.10) 



where G' is obtained from G by parallel (resp. series) expansion of an edge eo into a 
pair of edges ei, e2- It suffices to note that if |1 — p»| > 1 for z = 1, 2, then the same 
inequality holds for p\\ = 1 - (1 - pi)(l - p 2 ) and for Vseries = V1V2/ (pi + V2 - V1V2); 
the former is obvious, and the latter follows by observing that the series-expansion 
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formula corresponds to addition of 1/v = 1/p — 1 and that |1 — p\ > 1 corresponds 
to Re(l/f) < —1/2. In particular, since the generalized Brown-Colbourn conjecture 
manifestly holds for trees, it also holds for all connected graphs without a K± minor, 
as these are precisely the graphs that can be obtained from trees by a sequence 
of series and parallel expansions []37J, [72], |122| , 7JJ. The (original) Brown-Colbourn 



conjecture for series-parallel graphs was first proven by Wagner ||120|| , by a vastly 
more complicated method. 

4.2 Connected subgraphs containing a specified vertex 

Let G = (V, E) be a finite or countably infinite undirected graph equipped with 
edge weights {v e } ee E, and let x G V. Let us define the weighted sum over connected 
subgraphs G' = (V, E') C G containing n vertices, one of which is x, and m edges: 

C n>m (G,{Ve},x) J2 II ^ • ( 4 - H ) 

G' = (V',E') C G eG£' 
G' connected 
V 3 x 
\V'\ = n 
\E'\ = m 

Special cases are the tree sum 

T n (G,{v e },x) = C n ,n-i(G,R},x) (4.12) 
and the edge-counted sum 

m+1 

C., m (G,{v e },x) = J2 C n,UG,M,x) . (4.13) 

n=l 

When the edge weights v e are all equal to 1, we shall optionally omit them from the 
notation; note in particular the obvious bound 

C ntm (G,{v e },x) < C n>m (G,x) (sup\v e \) . (4.14) 

In this subsection we shall obtain a variety of upper bounds on C n , m (G, {v e },x) in 
terms of "local" information about the graph G and the weights {v e }. 

Proposition 4.2 Let G = (V, E) be a finite or countably infinite loopless undirected 
graph of maximum degree < r, equipped with edge weights {v e } eG E; and let x 6 V . 
Let T r be the infinite r -regular tree, and let y be any vertex in T r . Then 



C., m (G,x) < C., n {T r ,y) = T m+1 (T r ,y) = t% +1 = rLA^ (4.15) 



f(r) _ m [(r-l)(m + l)]! 

ml [(r — 2)m + r\\ 
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C. tm (G,{v e },x) < t% +1 (sup|u e |) . (4.16) 



ra-1 



and hence 



In particular, 

T n (G,{v e },x) < %> (sup\v e \ ] . (4.17) 

Proof. We can assume without loss of generality that G = (V, E) is connected. Let 
U = (V, E) be the universal covering graph of G, with covering map f:U^G; 
and let x be a vertex of U such that fix) = x. [The universal covering graph of a 
connected loopless graph G can be constructed as follows: Fix a base vertex x of G, 
and let the vertices of U be the walks in G (of finite length) that begin at x and do 
not contain any "doublebacks" (i.e. two consecutive uses of the same edge in opposite 
directions). Two vertices of U are defined to be adjacent if one of them is a one-step 
extension of the other, and f:U—>-G maps each walk onto its final vertex. We take 
x to be the zero-step walk starting at x.] It is easy to see that U is a tree (in general 
countably infinite even when G is finite). Moreover, since G has maximum degree < r, 
U is a subtree of T r , from which it follows trivially that C m>m (U,x) < C t>m (T r ,x). 
Let us prove, then, that C t>m (G,x) < C t>m (U,x). 

Fix an arbitrary total order on E, and choose arbitrarily for each edge e G E 
a distinguished direction. Now let if be a connected m-edge subgraph of G that 
contains x. Let S be the lexicographically first (with respect to the chosen total order 
on E) spanning tree of H. Then S based at x has a unique lifting to a subgraph S 
of U based at x: it is defined by mapping each vertex s of S to the unique path in S 
from x to s. 

Now, for each edge e of H not belonging to S, there is a unique edge e of U such 
that /(e) = e and e is incident with the image in S of the vertex of S from which e 
is directed. The addition of these edges to S produces a connected m-edge subgraph 
H of U that contains x. Moreover, the map H i— > H is injective, since H = f(H). 
This completes the proof that C m>m (G,x) < C, !m (T r ,x). 

We conclude by calculating the numbers t^ +1 = C, y . m (T r , x).f\ Let U r be the 
infinite tree in which all vertices have degree r except for one vertex y which has 
degree r — 1, and let u^ +1 = C, jm (U r ,y). Then define, as formal power series, the 
generating functions 

oo 

T r (z) = (4.18) 

n=l 

oo 

U r (z) = Y. U n )zn ( 419 ) 

71=1 



For similar computations, see e.g. p3[ 
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The recursive structure of r-regular rooted trees easily implies the functional equations 



T r (z) = z[l + U r (z)} r (4.20) 
U r {z) = z[l + U r {z)] r - X (4.21) 



We now use the Lagrange Implicit Function Theorem for formal power series ]48 



Theorem 1.2.4], which states that for formal power series f(u) = ^2^L Q f n u n and 
9( u ) = Yl^Lo 9nU n with g ^ 0, the functional equation U(z) = zg(U(z)) has a 
unique solution U(z), and for all n > 1 one has 

[z n ]f{U{z)) = l -lu n - l }{f{u)g{uY) (4.22) 
n 

where [z n ]P(2) denotes the coefficient of z n in the formal power series P(z). Applying 
this with f{u) = (1 + u) r and g{u) = (1 + u) r ~ l yields 

,(r) _ r f(r-l)n\ _ _ [(r-l)n]! 



tn n-l\ n-2 ) r (n- 1)1 [(r-2)n + 2]l ^' 2?>) 

(r) 1 f(r-l)n\ _ ^ [(r-l)n-l]! 



< J = - = (r - 1) t lv \ i F7 r J: — tt (4.24) 

n \ n-1 j v ; (n-l)![(r-2)n+l]l V ; 



Remarks. 1. The proof presented here is a simplification of my original proof, 
based on independent suggestions by Paul Seymour, Dave Wagner and an anonymous 
referee. 

2. A posteriori, we learn from Proposition ^]3|(c,d) below that the power series 
( |4.18| ) / ( |4.19| ) in fact define analytic functions in the disc \z\ < (r — 2) r ~ 2 /(r — l) r ~ J 



3. I suspect that Proposition fO| is known somewhere in the graph-theory liter- 
ature, but I do not know any reference. A weaker version of Proposition [4.2| can be 
found in [35, Lemma 5.4]. 



(r) 

Let us also collect some properties of the numbers tn that arise in Proposition 



4.2 



Proposition 4.3 The quantities 



t « = r [(r ~ 1)n]! (4 25) 

(n-l)![(r-2)n + 2]! ' 1 J 



defined for integers n, r > 1, have the following properties: 
_ / 1 for n = 1,2 



(a) n, 



for n > 3 
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(b) d 2) = n 

(c) As n — > oo at fixed r > 3, 
,(,) = r(r-l) 1 / 2 / (r-ir^ ' 

fd^ For a// n and a// r > 3 ; 



n 



-3/2 



1 _ _37_ _ 1 

r-1 r-2 



12n 



+ O 



t\ r) < f ^ K V ' < Kr-I)]- 1 



(r-2) 



r-2 



(4.26) 



(4.27) 



(e) As r — > oo at /ixed n > 1, 



(7j For a// r,n > 1, 



(rn 



,n-l 



' 3(ra- l)(n-2) 



2nr 



rn 



,71-1 



n! 



(4.28) 



(4.29) 



Proof, (a) and (b) are trivial, while (c) and (e) follow from Stirling's formula, (f) 
is trivial for r = 1, while for r > 2 it follows immediately from 



[(r-l)n]! 



(rn — n)\ 



< (rn) 



n-2 



[(r-2)n + 2]! (rn-2n + 2)! 
The first inequality in (d) is obvious for n = 1, so assume n > 2. We have 

rn 



(4.30) 



t 



(r) 



< 



< 



'(r — l)n 
[(r-2)n+l][(r-2)n + 2]A n 

r /(r — l)n 
(r — 2) 2 n \ n 

rfr ~ 1) 1/2 -a/a Ar-l) r - r 
27r(r- 2) 5 /2 \(r - 2) r - 2 

r(r-ir-4 „_ 3/2 /(r-l)«-^"- J 



2vr(r - 2) 



r+^ 



(r-2) 



r-2 



(4.31) 



where the second inequality uses Lemma O below. Then straightforward calculus 
shows that the function F(r) = r(r — l) r ~* /[y/2n(r — 2) r+ a] is decreasing on r > 
2; and we have F(3) = 12/v^F < 4 3 / 2 , F(4) = (27/8) < 3 3/2 and F(5) = 
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(1280/243) ^2/(371-) < 2 3 / 2 . So the first inequality in (d) follows except for the cases 
(r, n) = (3,2), (3,3) and (4,2), which can be checked by hand. The final inequality 
in (d) follows from 



lo g 7 7^37 = log a + (cr - 1) log 



(a-l)^ 1 ° v 7 a-l 



uu 



-A' 



' k(k + l) 



k=l 



= log a + 1 + log^l-^ (4.32) 
where the sums are convergent for a > 1, so that a "/ (er — 1) CT_1 < e(a — ~). ■ 



Lemma 4.4 Lei n >2 and 1 < k < n — 1 6e integers. Then 

n—k 



< 



n 



k) \n — k 



n 



2nk(n — k) 



(4.33) 



PROOF. We use the following strong form of Stirling's formula [j3T] , pp. 45-46]: for 
integer n > 1, 

logn! = (n + |)logn — n + logv / 2~7r + e n (4.34) 

with 

< e n < — . (4.35) 



12n 



12n + 1 

(The proof in [31] is valid only for n > 2, but ei = 1 — logv^vr ~ 0.08106 clearly 
satisfies 1/13 < e 1 < 1/12.) Then 



^k ^n—k *^ 



12n 12fc + l 12(n-Jfe) + l 

-144[n 2 - jfe(n - fc)] - 12n + 1 
12n(12fc + l) [12(n-fc) + l] 



< . 



(4.36) 



Proposition [12] clearly gives the best possible bound for C m>m (G,x) and T n (G,x) 
in terms of the maximum degree of G, since it is sharp when G = T r . On the other 
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hand, Proposition O is somewhat unnatural for general (unequal) edge weights {v e }, 
since adding an edge of small weight v e makes little change in C njm (G, {v e }, x) but can 
cause the bound to jump (in case it increases the maximum degree). It is of interest, 
therefore, to find alternative bounds that depend "smoothly" on the weights {v e }. 
We shall now give two such bounds (Propositions £|75| and |4~6| ). Unfortunately, both 
of them are strictly weaker than Proposition [O] when the edge weights are equal, 
and neither one is strictly stronger than the other. 

Proposition 4.5 Let G = (V, E) be a finite or countably infinite loopless undirected 
graph equipped with edge weights {f e } e gs- Then for any x G V , 

\ 

sup y^\v e \ (4.37a) 
^ v ^t ) 

(4.37b) 

[The e in front of the sup in (4 -37b) denotes, of course, the base of the natural 
logarithms!] 




Proof. As in the proof of Proposition |4.2|, we pass to the universal covering graph 



U = (V, E) of G with covering map f:U^G; and we define the weight v e of an edge 
e G E to be the weight v e of its image e = /(e). It then follows, as in Proposition 



0|, that C., m (G, {v e },x) < C. t JJJ, K},50- 

Let us now define, for each vertex x G V, the formal generating function 



C x (z) = ^C., m (f/,K},x)^ m . (4.38) 



m=0 

Then the recursive structure of rooted trees implies that 

C x {z) ■< X\[l + \v xy \zC y {z)} (4.39a) 



r< JJeK^M (4.39b) 



yr^x 

where y ~ x denotes that y is adjacent to x, xy denotes the (unique) corresponding 
edge, and ^ denotes coeflicientwise inequality at all orders in z\ the second inequality 
holds because 1 + az ^ e az for a > 0. It then follows, by induction on the power of 
z, that C x {z) -< C{z) for all x, where C(z) is determined by the equation 

C{z) = e» zC{z) (4.40) 
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with 

V" = sw P/.\ v xy\ = sw P/J v e\ ■ (4-41) 

tczV i&V ■ 

y^x e9« 

The coefficients of C(z) can be determined by applying the Lagrange Implicit Function 



Theorem to U(z) = zC(z), and we have the well-known (e.g. p. 392]) result 



m=0 



Finally, the inequality e m > (m + l) m /(m + 1)! is trivial for m = 0, 1, and for m > 2 
it follows from > x n /n\ by setting x = n = m + 1. I 

Remarks. 1. The proof presented here is a simplification of my original proof, based 
on the suggestions of an anonymous referee. 



2. Proposition [L5| holds also for graphs with loops, but they impose slight technical 
complications. 



An alternative estimate is due to Campanino et al. |29| , p. 129] (see also |28 



p. 522] and JTOl], pp. 463-464]): 



Proposition 4.6 Let G = (V, E) be a finite or countably infinite undirected graph 
equipped with edge weights {t> e } e6 £. Define the matrix M = [M xy ) x ^y by 

M xy = Yl H V2 • (4-43) 

e: e connects x to y 

Then for any x G V , 



2 m 



C. im (G,{v e },x) < (M 2m ) xx < sup^TM 1 / 2 . (4.44) 



iev 



Proof. Let G' = (V, E') be a connected subgraph of G having m edges. Then for 
any vertex x G V, there exists a path on G' starting and ending at x that uses each 
edge e G E' exactly twice. (Proof: The multigraph formed by doubling each edge of 
G' is Eulerian. Alternate proof: By induction on m|]) Conversely, every path on G 
starting and ending at x corresponds in this way to at most one subgraph G' . The 
claim follows. ■ 



Let us conclude by examining the relative sharpness of these bounds when G is 
an r-regular graph and the edge weights v e are equal. Then the tree bound tn ' of 

9 See e.g. fTM], Lemma V.7.A.2]. 
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Proposition [4.2| grows as n — > oo at an exponential rate (r — l) r ~V (r — 2) r ~ 2 : this is 
less than e(r — |) for all r, and behaves as e[r — | — 0(l/r)] as r — > oo. The bound of 
Proposition |4.5| is slightly weaker: it grows at exponential rate er. Finally, the bound 
of Proposition |4.6| grows at exponential rate r 2 , which is vastly weaker for large r but 
is slightly better when r = 2. 

In particular, when G is a regular lattice, it can be shown by supermultiplicativity 
arguments |BD|, ^T], |127| , p7| ] that the limits 

X (G) = lim T n {G,x) 1/n (4.45) 

n— +00 

X b {G) = lim C., m (G,x) 1/m (4.46) 

m— ->oo 

exist. For the simple hypercubic lattice 7L d with nearest-neighbor bonds, these growth 
constants have been computed (non-rigorously) in a large-d asymptotic expansion 
J5L U7\ (see also |7S|, |5S], for related rigorous results): 

logA G (Z d ) = log a + 1 - V 1 - *-a~ 2 - ... (4.47) 

log A fc (Z d ) = log a + 1 - ^a- 1 - f ~ - a" 2 - . . . (4.48) 

where a = r — 1 = 2d — 1. Let us compare this with the tree bound of Proposition 
L|: 

Thus, the latter bound is very close to sharp for G = Z d in high dimension d, con- 
firming the intuition that high-dimensional regular lattices are "like trees" to leading 
order in 1/d. 



5 Application to the Potts-Model Partition Func- 
tion 

We are now ready for the main theorem of this paper: 

Theorem 5.1 Let G = (V,E) be a loopless finite undirected graph equipped with 
complex edge weights {v e } e£ E satisfying \l+v e \ < 1 for all e. Let Q = Q(G, {v e }) > 
be the smallest number for which 

00 

inf a' 1 "pe an Q-( n -^ m&xT n (G,{v e },x) < 1. (5.1) 

a>0 ^ x£V 
n=2 

[Note that Q is automatically finite, since T n (G, {v e },x) = for n > \V\.J Then all 
the zeros of Z G (q, {v e }) lie in the disc \q\ < Q. 



24 



Proof. Starting from the polymer-gas representation (|2.2|) - (l2~4]) of Zc(q, {^e}), we 
apply Theorem |3.1| and Proposition |3.2| with the choice R$ = \w(S)\. We verify 
hypothesis (a) of Proposition [T2] by verifying ( |3.15| ) with 

An = max V \w(S)\ . (5.2) 

\S\=n 



Now we use Proposition ^4] to conclude that w(S) can be bounded by a sum over 

\w(S)\ < \q\-^ Yi nw- ^ 

EC E ee-B 
(S,E) tree 



Inserting this into (|5.2j) , we get 

A n < |gr (n - 1} maxT n (G,{v e },x) . (5.4) 



If |g| > Q, hypothesis (b) of Proposition |3.2| holds (recall Remark 4 following that 
Proposition) and hence Zc(q, {v e }) 7^ 0. I 



In applying Theorem |5.1| we are of course free to use any convenient upper bound 
on ma.x x( zv T n (G , {v e },x). In particular, when G has maximum degree < r, Proposi- 
tion [42| provides such a bound. Recall that 

t« = r [{r - 1)n]l (5 5) 

n " (n-l)![(r-2)n + 2]! ' K } 

and let C = C(r) > be the smallest number for which 

oo 

inf a' 1 Ve an C- (n -^4 r) < 1 . (5.6) 

n=2 



The following is then an immediate consequence of Theorem 5.1 and Proposition 4.2 



Corollary 5.2 Let G = (V, E) be a loopless finite undirected graph of maximum 
degree < r, equipped with complex edge weights {v e } e€ E satisfying \l + v e \ < 1 for alle. 
Let v max = max \v e \. Then all the zeros of Zc{q, {v e }) lie in the disc \q\ < C(r)v max . 

e£E 

And for the chromatic polynomials: 

Corollary 5.3 Let G = (V, E) be a loopless finite undirected graph of maximum 
degree < r. Then all the zeros of Pa{q) tie i n the disc \q\ < C(r). 
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Table |l] lists rigorous upper bounds on C(r) for 2 < r < 20, proven (with the as- 
sistance of Mathematica) as follows: After computing numerically an approximate 
value of C(r),F] I added 10~ 6 and rounded it upwards to a rational number p/10 6 . 
(Thus, the value reported in Table p] exceeds my best estimate of C{r) by at most 
2 x 10 -6 .) I likewise approximated the numerically-found a by a rational number 
p'/10 6 . Thereafter I did all computations in exact rational arithmetic. First I com- 
puted a rational upper bound on e a (differing from the true e a by at most 2 x 10~ 10 ) 
by truncating the Taylor series for e~ a at odd order (here ninth or eleventh) to obtain 
a lower bound on e~ a . Finally, I computed an upper bound on ( |5.6|) by summing the 
terms explicitly through n = some no and bounding the tail of the series (n > n$ + 1) 



using Proposition fO](d); I systematically increased n until the inequality ( |5.6|) was 
verified. For r = 2, of course, I just summed the series exactly. 
As r — > oo we have the following: 

Proposition 5.4 Let K ph 7.963906. . . be the smallest number for which 

oo n _i 

inf a' 1 e an K- {n - 1] ^— < 1. (5.7) 



n\ 

n=2 



Then C(r) < Kr for all r, and lim C(r)/r = K. Moreover, we have the rigorous 

r— >oo 

bound K < 7.963907. 



PROOF. Clearly C(r) = C(r)/r is the smallest number for which 

oo ^ Jr) 

inf a" 1 Y e an ^- < 1 . (5i 

n=2 



It follows from Proposition [4.3| (f) that C(r) < K for all r. 

Now suppose that we were to have lim inf C(r) < K — e < K. Then there would 

r— >oo 

exist infinite sequences {r{\ j oo and {a{\ such that 



^e^iK-e)-^^ < 1 (5.9) 



* — . r 

n=2 



for all i. Now the finiteness of ( |5.9| ) implies that the oti are bounded [e.g. from 
Proposition |4.3|(c) we have e ai < (K — e)r i (r i — 2) ri ~ 2 / (ri — iy^ 1 < j(K— e) whenever 

10 Using the generating function f(z) — Y^=2^' >z ™ wnere z = e^/Cj I solve simultane- 
ously the equations f(z) = (logz + logC)/C7 and f'(z) = l/(Cz) by solving numerically f(z) = 
—zf'(z) log/'(z) and then plugging back in to determine C = l/[zf'(z)} and a = — log/'(z). 
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20 
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28 



Table 1: Upper bounds on C{r) for 2 < r < 20; they differ from my best estimate 
of the true C(r) by at most 2 x 10 -6 . The third column gives a value of a for which 
( |5.6|) is proven to be < 1. The fourth column gives the upper bound on e a employed 
in this proof. The fifth column gives the number of terms explicitly summed in the 
series. 
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r,i > 3]. So we can extract a subsequence of that converges to some value a*. 
Then Proposition |4.3| (e,f ) and the dominated convergence theorem imply that 



OO „_1 

n n 

a 

n 

n=2 



n 

Te^fif-eJ-^'V < 1, (5-10) 



which contradicts the definition of K. 

Finally, it is easy to prove that K < 7.963907, by a computer-assisted method 
similar to that used above for C(r). For the tail of the series {n > no + 1), it suffices 
to use the crude bound n n_1 jn\ < e n ~ l < 3 n . The proof succeeds with the choices 
a = 0.403774, e a < 1.4974655 and n = 32. I 



If we employ Proposition [4.5| in place of Proposition [0|, Theorem |5.1| yields the 
following: 

Corollary 5.5 Let G = (V, E) be a loopless finite undirected graph equipped with 
complex edge weights {v e } e& E satisfying |1 + v e \ < 1 for all e. Then all the zeros of 
Za(q, {v e }) He in the disc \q\ < K maxj^ \v e \, where K is the constant defined in 

Proposition |5-4 

Remarks. 1. Since Zc(q, {v e })/q for any graph G is the product of the same 
quantity over the blocks of G, it is legitimate to apply Theorem |5.1| and its corollaries 
separately to each block. This can lead to large improvements (consider e.g. trees). 

2. What happens if we drop the assumption that |1 +v e \ < 1? Because we can no 
longer use Proposition |4.1| to reduce the sum to trees, we need to consider all n- vertex 
connected subgraphs of G containing a given vertex x. But the number m of edges 
in such a subgraph could be as large as [rn/2\ (where r is the maximum degree of 
G). Therefore, the factor tn^v^l coming from Q4.171 ) has to be replaced by a factor 
^m=n-i^m+i v max coming from (|4.16| ) [or the analogue from Proposition |4.5| . As a 



consequence, the radius of the g-plane disc containing all the zeros of Zc(q, {v e }) 
will scale as m&x[v max , Vmax] rather than simply v max . And this is not simply an 
artifact of the method of proof: for the g-state Potts ferromagnet (v > 0, q > 0) 
on the the simple hypercubic lattice 7L d with nearest-neighbor bonds, the first-order 
phase-transition point v t indeed behaves as 

v t (q) = ^[1 + 0(1/9)] (5.11) 

as q — > +oo |73, |67|, |64|, |66|, |16|]. In the Whitney-Tutte-Fortuin-Kasteleyn polynomial 



( |1.2|) , this reflects the coexistence at v = v t (for all g > 1) between a phase with a 
low density of occupied edges and a phase with a high density of occupied edges. 
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6 Some Generalizations 

The following generalization of the Whitney-Tutte-Fortuin-Kasteleyn polynomial 
(|0]) is motivated by some work of Tutte J[T§ JTT|], Farrell |KJ and Stanley |IU§ |TD§ 
as well as by the statistical-mechanical application to be discussed below. Let us 
replace the single complex number q by a map q: — > C, and define 



Z G (q,{ve}) = £ ( nq(V0] 



(6.1) 



where (Vi, Ei), . . . , ^K e ')) are ^ ne connec t e d components of (V, E'). We im- 



mediately deduce an analogue of Proposition |2.1| : the identity ( |2.2| ) is replaced by 



(6.2) 



and the fugacities w(S) are now given by 

q(g) 

u;(5) 1 



n , , E if|5|> 2 

Hq(M) ECE e&E 

x£S (S, E) connected 



(6.3) 



if |5| < 1 

The proof of Theorem |5.1| then goes through without change, and yields: 



Theorem 6.1 Let G = {V,E) be a loopless finite undirected graph equipped with 
complex edge weights {v e } e€ E satisfying |1 + v e \ < 1 for all e, and let q: V*(V) — > C. 
Let {Q n }^ =l be a sequence of positive numbers satisfying 



inf a^r^Q: 1 maxTJG, {v e },x) < 1 

Ct>0 ' J itCV 



n=2 



and assume that 



q(S) 



n q({^}) 



— Qlsj 



(6.4) 



(6.5) 



for all nonempty subsets S C V . Then Zg(c\, {v e }) 7^ 0. 

The following special case is of particular interest: Fix an integer N > 0, and for 



each x G V choose a vector u 



- K u x u x N ^) e C w . Then define 

A/ 



q(5) = q - N + J2l[(l + 



u 



(^)^ 
a; / 



(6.6) 
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where q is a fixed complex number. This corresponds to a g-state Potts model in 
a magnetic field h x = ( hi 1 hx) in the first N spin directions, where u$ = 
exp(h x l) ) - 1. To see this, we first define, for each integer q > N, the partition 



function for the 9- state Potts model in a magnetic field, generalizing (O): 



Zg(q> R}> {4°}) 



{a x } e&E 



x&V i=l 



Now expand out the product over e G £?, and let £" C be the set of edges for which 
the term taken. Then perform the sum over configurations {<J X }: in 

each connected component of the subgraph (V, E') the spin value a x must be constant, 
and there are no other constraints. The sum over possible spin values in a connected 
component with vertex set S yields ( |6.6| ). It follows that, for any integer q > N, 
the partition function Zc(q, {v e }, {u x }) equals Zc(q, {v e }) with weights (|6.6|). We 
then take the latter, which is a polynomial in q, {v e } and {u x ^}, as the definition of 
Zo(q, {v e }, for general complex q. 

The following lemma gives a sufficient condition for the applicability of Theorem 



5.1 to this situation: 



Lemma 6.2 Let q(S) be defined by ([^7 
(a) If N = 1 and each u* x equals either or —1, then 



q(s) 



EI q(W) 



< min(|g|, |g - 1|)- (|5| ~ 1} 



(6.8) 



(b) If — 1 < u x ^ < for all x,i and \q\ > N, then 



q(S) 



n q(W) 

xGS 



< Q q \_ N y(\s\-D 



(6.9) 



(c) If |1 + ui | < 1 for all x, i and \q — N\ > N, then 



q(s) 



n q({^» 

xes 



< 



\q-N\+N 
(\q — N\ — N)\ s \ 



(6.10) 



The proof of Lemma |6.2| is deferred to the end of this section. 
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We can exploit this example to obtain new results for the ordinary (zero-field) 
Potts-model partition function Z G (q, {v e }) and in particular for the chromatic poly- 
nomial P G {q), by employing a variant of the "ghost spin" trick of Suzuki [|110|| and 



Griffiths j49|. Given a finite graph G = (Vo,E ) an d an integer N > 1, we de- 
fine G to be the join of Go with the complete graph on iV vertices. Thus, the 
vertex set of G is V = Vq U{2/i; • • • j Vn} (disjoint union) and the edge set is E — 
^oU{( :E ^)}xGVb,i<i<7vU{(?/i%')}i<i<j<Af- We allow the ed S e weights {v e } eeEo and 
{ v (xyi)}xev , i<i<N to be arbitrary complex numbers, but we require that = — 1 

for 1 < i < j < N (this condition is crucial). We then have the identity 

Z G (q,{v e ,v {xyi) ,v {y . yj) }) = q {N) Z Go (q,{v e },{u^}) , (6.11) 

where q/m = q(q — 1) • • • (q — N + 1) is the iV-th "falling factorial" polynomial and 

ui 3 = Vtxyi)- This is most easily proven in the Potts spin representation (|1.1| )/ (|6\7| ): 
Let q be an integer > N, and let us compute the left-hand side of (|6.11|) . There are 



g(7v> admissible ways to color the vertices {yi, . . . ,Dn}, all of which are equivalent 
modulo permutations of {1, . . . , q}; and with any such coloring fixed, the sum over 
colorings of V yields precisely Z Go (q, {v e }, {u^}) with Ux = v {x Vl )- Since both sides 
of (|6.11| ) are polynomials in q and the equality holds for infinitely many values of q, 



it must hold identically. 

By applying Theorem |6.1| to the graph Go, we can obtain new results for the 
ordinary Potts-model partition function of the graph G. In particular, given any 
graph G = (V, E) and any vertex y 6 V, we can interpret G as the join of G = G \ y 
(the graph obtained from G by deleting y and all edges incident on it) and K\. (Any 
edge {xy) that was not originally present in G can be introduced and given v^ xy ) = 0.) 
More generally, given any iV-clique yi, . . . , of G, we can interpret G as the join of 
G = G \ {y%, . . . , yjy} and K N ; however, for iV > 1 we must require that v ( yiy .) = — 1 
for each pair i ^ j. Theorem |6.1| , Lemma |6.2j (b,c) and Proposition [4.2| then yield 



an extension of Corollary |5.2| . To state it, we first define G = C(r,N,v) to be the 
smallest number for which 



00 c _i_ /v 

inf a' 1 V e an 1)"" 1 ^ < 1. (6.12) 

a>0 ^ {C _ N) n ' - 



We then have: 



Theorem 6.3 Let G = (V, E) be a loopless finite undirected graph in which all ver- 
tices have degree < r except perhaps for an N -clique yi, . . . ,yN- Let G be equipped 
with complex edge weights {v e } e€ E satisfying |1 + v e \ < 1 for all e and V( y . y .) = — 1 
for all i ^ j . Let v max = max \v e \, where E is the set of edges not incident on any of 

the vertices y%, . . . ,y^- Then: 
(a) All the zeros of Z G (q, {v e }) He in the disc \q — N\ < G(r, N,v max ) . 
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(b) If, in addition, all the edges e incident on any of the vertices yi, . . . ,Un satisfy 
— 1 < v e < 0, then all the zeros of Z G (q, {v e }) lie in the disc \q\ < C(r)v max +N. 

And for the chromatic polynomials, we have, using Lemma |B~2] (a): 

Corollary 6.4 Let G = (V, E) be a loopless finite undirected graph in which all 
vertices, except perhaps one, have degree < r. Then all the zeros of Pg(q) tie i n the 
union of the discs \q\ < C(r) and \q—l\ < C(r). In particular, they all lie in the disc 
\q\ < C(r) + 1. 

Thus, the zeros of Pg(q) can be bounded in terms of the second-largest degree of 



a vertex in G. Such a result was recently conjectured by Shrock and Tsai |p9| ; see 
Section [7| for further discussion. 

Let us note that the phrase "except perhaps one" in Corollary |6]4] cannot be 
replaced here by "except perhaps two", not even in the case r = 2. Indeed, I have 



elsewhere ||103|| constructed a family of planar graphs in which all but two vertices 



have degree 2 and whose chromatic roots are together dense in {q G C: \q — 1| > 1}. 



Modifications of these graphs show also ||103|| that the condition v^ yiyj ) = —1 for i ^ j 



in Theorem |6.3| (when N > 1) cannot be relaxed. 



Let us now give the proof of Lemma |6.2| . We will need the following elementary 
fact: 

Lemma 6.5 Let z and a be complex numbers. Then 

\z + Xa\ 2 > \z + a\ (\z\ - \a\) (6.13) 

whenever < A < 1. 

Proof. Simple calculus shows that 



r|.:| 2 ifRe(^*a)>0 

* „\2 



min \z + Aal 

0<A<1 



A 2 - Re | y if-|a| 2 < Re(z*a) < (6.14) 
a 



Jz + a| 2 if Re(z*a) < -\a\ 2 

In the first two cases we clearly have 

min \z + Xa\ 2 > \z\ 2 - \a\ 2 = (\z\ + \a\)(\z\ - \a\) > \z + a\ (\z\ - \a\) , (6.15) 

0^ 1 

while in the third case we have 

min b + Aal 2 = \z + a| 2 > \z + a\ (\z\ — \a\) . (6.16) 

0<A<1 
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Proof of Lemma We use the shorthand it? (5) = q(S)/ ]J q({x}). 

xes 

(a) Let \S\ = n and suppose that the sequence (ui 1 )xes consists of m — l's and 
n — m O's. Then 



w(S) 



Q 



-(n-1) 



if m = 



(6.17) 



from which ( p.S\ ) immediately follows. 

(b) Let S = {xi, . . . , £„}; we then have q({xj}) = q + Uj and q(S) = q + u with 
— < u < Ui, . . . , u n < 0. Now apply Lemma |0| with z = q, a = u and A = 
for j = 1, 2: we have 



q(5) 



q({^i})q({^2}) 



i i 
< — — — < 



|g| — |it| 



|g| — iV 



and hence 



1 1 
g - N - L -^ 5 + it, 



(6.18) 



(6.19) 



which implies ( p.9| ). 

(c) This bound is trivially obtained by bounding the numerator and denominator 
separately. I 



Remarks. 1. For simplicity, I have not bothered to exploit the full strength of Q6.19Q , 
which is quite a bit sharper than (|6.9[) . 



2. I am not entirely happy with Lemma p.2| , and I suspect that it can be improved. 
In particular, it is disconcerting that (|6.9| ) is not uniformly stronger than (|6.10|) , even 



though the corresponding hypothesis on the Ux is strictly stronger 



7 Some Conjectures and Open Questions 

The bounds in this paper are, of course, far from sharp, and it is of some interest 
to speculate on what the best-possible results might be. Let us define 



Cnnt(r) 



j o P t\i ) — max{|g|: Pg(q) = for some loopless graph G of maximum degree r} . 



(7.1) 
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The example of the complete graph K r+ i shows that C opt (r) > r. It is easy to see that 
C ptiX) = 1 and C opt {2) = 2; and there is some evidence that C opt (3) = 3.0 But, at 
least for r > 4, C opt (r) must in fact be strictly larger than r, as is shown by numerical 
computations on the complete bipartite graph K r r (see Table Indeed, Gordon 
Royle (private communication) has conjectured that, for r > 4, K r ^ r is the graph 
of maximum degree r having the largest chromatic roots (in modulus). It would be 
useful to have a better understanding of the chromatic zeros of the complete bipartite 
graphs K m ^ n . In particular, it would be useful to have a proof that K r r has chromatic 
roots of magnitude > r for all r > 4; and it would be valuable to understand the 
asymptotic behavior of the chromatic roots of K m>n as m, n — > oo in various ways 
(e.g. with a = m/n fixed). 

Using the Dobrushin uniqueness theorem p6|, |101|| , it can be proven p9| that for a 



countable graph G of maximum degree r, the g-state Potts-model Gibbs measure on 
G is unique for all integer q > 2r whenever — 1 < v e < for all edges e. Uniqueness 
of the Gibbs measure is one of several (inequivalent) notions of "absence of phase 
transition" |101| . It does not imply the analyticity of the free energy, but it does 
make it plausible.^ Likewise, a result that holds for integer q > qo need not hold for 
all real q > qo, much less for a complex neighborhood of that real semi-axis; but it 
does suggest that such a result might be true. It is not unreasonable, therefore, to 
conjecture that there is a complex domain D r containing the interval (2r, oo) of the 
real axis, such that Zo(q, {v e }) ^ whenever q G D r , — 1 < v e < for all edges e, 
and G has maximum degree < r. Indeed, it is quite possible that D r = {q: \q\ > 2r} 
works; this would be a slight extension of the conjecture that C opt (r) < 2r. 

We can pose these questions more generally as follows: Let Q be a class of finite 
graphs, and let V be a subset of the complex plane. Then we can ask about the sets 

Si(G,V) = |J |J {qeC: Z G (q,v) = 0} (7.2) 
Gee vev 

S*{G,V) = |J |J {q G C: Z G (q, {v e }) = 0} (7.3) 
Geg {v e }-. tj e ev Ve 



11 Biggs, Damerell and Sands |l2| have verified that the chromatic roots of all 3-regular graphs 
with < 10 vertices, as well as those of ladders ("prisms") and Mobius ladders of arbitrary length, 
lie in \q\ < 3. Read and Royle |8^] have extended this verification to all 3-regular graphs with < 16 
vertices, as well as to some larger graphs. 



12 Recall (m|, g that 

m 



k=0 

where S(m, k) is the Stirling number of t he se cond kind (the number of ways of partitioning a set 



of to ele ment s into k nonempty subset s) [105 , pp. 33-38] and q^) = q(q — 1) • • • (<jf — + 1). See 
Woodall |l28| , pp. 219-220] and Brown |l8| for some properties of the chromatic zeros of the K m . n . 

13 Indeed, it was by meditating on possible extensions of the theorem in [^9| that I was led to 
conjecture the results in this paper. 
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r 



q 



2 


1.500000 ± 0. 


3 


2.140640 ± 1. 


4 


2.802489 ± 3. 


5 


3.469365 ± 4. 


6 


4.138450 ± 5. 


7 


4.808805 ± 6. 


8 


5.480007 ± 8. 


9 


6.151830 ± 9. 


10 


6.824136 ± 10. 


11 


7.496833 ± 11. 


12 


8.169855 ± 13. 


13 


8.843156 ± 14. 


14 


9.516697 ± 15. 


15 


10.190450 ± 17. 


16 


10.864391 ± 18. 


17 


11.538501 ± 19. 


18 


12.212764 ±21. 


19 


12.887165 ±22. 


20 


13.561693 ±23. 



866025 % 1.732051 

948682 % 2.894772 

097444 % 4.177093 

291184i 5.518221 

516667 % 6.896404 

765768 i 8.300616 

033190 i 9.724331 

315289i 11.163316 

609446 i 12.614641 

91371H 14.076186 

226591 i 15.546358 

546915 % 17.023928 

873744 % 18.507925 

206318 i 19.997566 

544006 i 21.492211 

886280 % 22.991328 

232697 i 24.494469 

582876 % 26.001256 

936489 % 27.511362 



Table 2: The chromatic root of largest modulus for the complete bipartite graphs 
K r , r for 2 < r < 20. 
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Among the interesting cases are the chromatic polynomials V = { — 1}, the antiferro- 
magnetic Potts models V = [—1,0], and the complex antiferromagnetic Potts models 
V = A = {f £ C: < 1}. Indeed, one moral of this paper is that some questions 

concerning chromatic polynomials are most naturally studied in the more general 
context of antiferromagnetic or complex antiferromagnetic Potts models (with not- 
necessarily-equal edge weights). In Corollary |572| we have shown that the set S2(Q r , A) 
is bounded, where Q r is the set of all loopless graphs of maximum degree < r; and 
in Theorem |T3] we have extended this to S 2 (Q' r , A), where Q' r is the set of all loopless 
graphs of second-largest degree < r. But it would be interesting to examine in more 
detail the location of all these sets in the complex plane, and to prove sharper bounds. 

Another direction in which the results of this paper could be extended is by finding 
a criterion weaker than bounded maximum degree (or bounded second-largest degree) 
under which the zeros of Pg(q) an d Z G (q, {v e }) could be shown to be bounded. An 
interesting idea was suggested very recently by Shrock and Tsai pSj] , who studied 
a variety of families of graphs and arrived at a conjecture that can be rephrased as 
follows: For G = (V, E) and x, y e V, define 

X(x,y) = max # of edge-disjoint paths from x to y (7.4a) 

= min # of edges separating x from y (7-4b) 

and 

A(G) = maxX(x,y) . (7.5) 

Clearly X(x,y) < min[deg(x), deg(y)] and hence A(G) < second-largest degree of G. 
Now let Q A be the set of all loopless graphs with A(G) < r. Then the conjecture is 
that the set S 2 {G A ,V) is bounded, where V = { — 1} or [—1,0] or perhaps even 
More generally, one could define X(x, y; {v e }) to be the maximum flow from x to y 
when \v e \ is taken to be the capacity of edge e, and likewise A(G, {v e }); this might lead 
to the appropriate extension of Corollary |5.5| . This possible connection of chromatic- 
polynomial and Potts-model problems with max-flow problems is intriguing. Note 
that A(G) and A(G, {v e }) possess a "naturalness" property that maximum degree 
and its relatives lack: namely, for any graph G with blocks G\, . . . ,Gb, we have 
A(G, {v e }) = max A(Gj, {v e }); contrast this with Remark 1 after Corollary |5.5| . 

Ki<b 



14 Shrock and Tsai |9^| studied only the chromatic-polynomial case V = { — 1}, and proposed an 
even stronger result, based on the quantity 

Kon-adj(G) = max X(x,y) . 

x, y not adjacent 



But this cannot work for u ^ -1: a counterexample is obtained |103| by gluing together n copies of 
the cycle Ck (any fixed k > 3) along a single common edge and then taking n — > oo. 
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